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Abstract 

A mathematical formalism for the analysis of correlations in multi-source events 
such as W + W~ production in e + e~ annihilations is presented. Various measures 
used in experimental searches for inter-VF correlations are reviewed. 



1 Introduction 



The problem of possible inter-VF dynamical correlations continues to attract much 
experimental attention at LEP-II []J . The importance of these in a precision measure- 
ment of the H^-mass ||-ff] is by far the main reason for the flurry of present activities. 
However, this should not be the sole motivation for careful experimental work in this 
subfield. For example, the possible absence of Bose-Einstein correlations (BEC) be- 
tween pions originating from different W's raises interesting general and basic questions 
regarding the coherent versus incoherent nature of particle emission. Also the question 
of possible color- or string-reconnection effects may be of considerable importance for 
our understanding of the vacuum properties of Quantumchromodynamics. 

As repeatedly emphasized by the Lund group ]5|,|[, Bose-Einstein correlations of 
a coherent type, for which we suggest the name "String Symmetrization Correlations 
(SSC)", are present in any string model of hadron production. Moreover, the SSC 
depend essentially only on local properties of the string and should thus be indepen- 
dent of the environment in which the string is fragmenting. This basic property is 
not in contradiction with experimental results on Bose-Einstein correlations in e + e~ 
annihilations and lepton-nucleon scattering. 

However, in systems comprising several strings, a second type of Bose-Einstein cor- 
relations may exist if the different strings behave independently so that they act as 
incoherent sources of particle emission. This second-order intensity correlation effect, 
or HBT effect, since first discovered by Hanbury Brown and Twiss 0, is expected 
to reflect the "geometry" of the collision process and, in particular, to be sensitive 
to the size of the "freeze-out" volume where the hadrons are formed. This volume 
could have an extension of many fermi's and can thus be significantly larger than the 
typical "radii" of less than a fermi, commonly measured e.g. in BEC studies of e + e~ 
annihilations. As a result, correlation functions of identical bosons may well show an 
additional enhancement at substantially smaller values of their momentum difference 
than for SSC correlations. 

The reaction e + e~ — > W + W _ — > qiq2q3q4 — > hadrons is a prototypical example of 
a system of which the hadronisation proceeds via the fragmentation of two color-fields 
or strings. Barring possible reconnection effects at an early stage of the evolution of 
the system, these strings are thought to fragment independently. Thus, SSC within 
single W's and inter-VT HBT correlations could coexist in this system, but manifest 
themselves in different ranges of e.g. the commonly used variable Q 2 , the square of 
the difference of the four-momenta of the identical pions. No dedicated searches in 
this direction have been made so far. Moreover, if the HBT correlations are of much 
shorter range than naively expected, the influence on IV-mass measurements is likely 
to be much weaker than present studies suggest. 

Limited statistics as well as limited experimental sensitivity at very small Q may 
prevent the HBT-type of correlation effect to be clearly observed in e + e~ — > W + W~ — > 
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Aq events. However, as yet unexplored alternative reactions exist. Much larger statistics 
is available in e.g. e + e~ —>■ qq + g — > three-jet events. For these, due to the emission of 
a hard gluon, two color-disconnected strings are stretched, one between the quark and 
the gluon, another between the gluon and the antiquark. Depending on the origin of 
the identical pions studied, SSC as well as HBT correlations should contribute^] to the 
correlations among identical bosons. 

In hadron-hadron, hadron-nucleus and nucleus-nucleus collisions, hadroproduction 
is also believed to result from the break-up of several to very many strings. Since the su- 
perposition of many independent particle "sources" weakens the measured strength of 
correlations of the SSC type, HBT correlations may well dominate the Bose-Einstein 
correlations measured in these processes. This could also explain the observed pos- 
itive correlation between particle density (or multiplicity) and the measured BEC 
radii. Finally, due to the "inside-outside" character of hadron production, whereby 
low-momentum particles "freeze out" first (in any reference frame), one may expect a 
correlation between the width of the Bose-Einstein enhancement and particle momenta 
if the HBT correlations are the dominant effect. 

Returning to correlations in the WW system, experimental study of such effects 
should start from a well-defined mathematical framework. Since, in general, inter-VT 
dynamics introduces genuine correlations between the decay products of the W's, we 
generalize the formalism presented in an earlier paper on the subject |§ where observ- 
ables relevant for the case of stochastic independence and fully overlapping decays of 
the Ws were proposed. We consider the general case of stochastic dependence and the 
separation of the W hadronic decay products in momentum space. We concentrate on 
second-order (two-particle) inclusive densities and correlation functions but the results 
can be generalized to higher orders. 

In the mathematical treatment of the problem, a general inter-py correlation func- 
tion is introduced, representing an arbitrary stochastic correlation between hadrons 
from different Ws which could arise from color reconnection effects, Bose-Einstein 
correlations or others. Although most of the examples treated in this paper relate to 
Bose-Einstein studies involving identical particles, the formalism is general and can be 
used, where needed with a suitable change of kinematic variables, in other contexts as 
well. 

2 Multivariate distributions and moments 

Before discussing the problem of correlations among particles originating from W 7 s 
decaying into fully hadronic, so-called four-jet configurations, e + e~ — > qiq 2 q3Q4 — ► 
hadrons, we first consider a more general problem. 

Assume that a system (e.g. an event) comprising in total n (observed) particles, 

1 I thank B. Andersson and G. Gustafson for discussions of this point. 
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can be subdivided into S possibly stochastically correlated groups or "sources" f2j 
(i = 1, . . . , S). We take the n particles to be of identical type and assume that the 
groups are mutually exclusive i.e. every particle can be assigned to one and one group 
onlyQ: Vti nQj = 0, and i ^ j. No assumptions are made, however, about possible 
overlap in momentum space of particles from different groups. 

Let n be the total number of particles counted in the union of the 5* groups, 



s 

n= J2 n m ■ (1) 

m=l 

Consider further the multivariate multiplicity distribution P$(ni, . . . , ns) giving the 
joint probability for the simultaneous occurrence of rii particles in class Qj., . . . , n$ 
particles in class Q$- The probability distribution of n is then given by 

n n 

p ( n ) = Z) ' " £ p s{ni, n s )5 ntni+ ... +ns . (2) 

ni=0 ng=0 

We further define the single- variate factorial moment generating function of P(n) 

oo 

G{z) = "£{l + zrP{n). (3) 

n=0 

By expanding (1 + z) n , (|3|) can be rewritten as 

gw = ^Sfii^V^ (4) 

g(z) = i + E^r^ (6) 

g=i y- 

where in the last line we introduce the symbol F q for the (unnormalized) factorial (or 
binomial) moment of the probability distribution P(n). The brackets in Eq. (^) denote 
a statistical average. 

F q = (n(n-l)...(n-q + l)), 

= / dj/i... / dy q p q (y u . . . ,y q ). (7) 

J A J A 



2 Such a subdivision can be based on a "natural" partition of the n particles in S groups according 
to the underlying dynamics. It could also be based on an experimentally dictated partition of the 
phase space in S distinct (non-overlapping) regions e.g. as the result of jet clustering. 



3 



The last equation expresses that, for identical particles, the factorial moment F q is 
equal to the integral over g-dimensional phase space (here for simplicity of notation 
represented by the variables yi) of the g-particle inclusive density p q (yi, ■ ■ ■ ,y q ) over 
the same phase space volume A ||. The g-particle inclusive density p q (yi, . . . ,y q ) is 
defined as 

, . 1 d q a ... 

^•••'^ ) = ^' (8) 

with a the total cross section of the considered reaction. Experimentally, this quantity 
is approximated by 

1 diV^P 168 

with fijV^uples the number of g-tuples of particles, counted in a phase space domain 
(yi + dyi, . . . , y q + dy q ); N ev ^ is the number of events in the sample. 

Factorial cumulants are formally defined as the coefficients of z q /q\ in the Taylor 
expansion of the function \ogG(z): 

°° z q ~ 

\ogG(z) = (n)z + Y / ^K q . (10) 

q=2 1- 

The (unnormalized) factorial cumulants of order q, K q , also known as Mueller 
moments |TUJ, are equal to the g-fold phase space integral of the g-particle inclusive (so- 
called "connected" or "genuine") factorial cumulant correlation function C q (y\, . . . , y q ) 

^ q = I dyi " ' Ia dyq Cq ^ yu ( n ) 

The correlation functions C q (yx, . . . , y q ) are, as in the cluster expansion in statistical 
mechanics, defined via the sequence 

Pi(l) = 0,(1), 
p 2 (l,2) = C 1 (1)C 1 (2) + C 2 (1,2), 
p 3 (l,2,3) = C 1 (1)C 1 (2)C 1 (3) + C 1 (1)C 2 (2,3) + C 1 (2)C 2 (1,3) + 

+ C 1 (3)C 2 (1,2)+C 3 (1,2,3), (12) 
etc. (13) 

These relations can be inverted yielding 

C 2 (l,2) = p 2 (l,2)- Pl (l) Pl (2) , 
C 3 (l,2,3) = p 3 (l,2,3)-^p 1 (l)p 2 (2,3) + 2p 1 (l)p 1 (2)p 1 (3) , 

(3) 

C 4 (l,2,3,4) = p 4 (l,2,3,4)-^ Pl (l)p 3 (l,2,3)-^p 2 (l,2)p 2 (3,4) 

(4) (3) 
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+ 2^:p 1 (l)p 1 (2)p 2 (3,4)-6p 1 (l)p 1 (2)p 1 (3)p 1 (4), 

(6) 



etc. 



(14) 
(15) 



In the above relations we have abbreviated p q (yi, ■ ■ ■ , y g ) to p g (l, 2, . . . , q) etc.; the 
summations indicate that all possible permutations have to be taken (the number 
under the summation sign indicates the number of terms). 

The factorial moments F q and factorial cumulants K q are easily found if G(z) is 
known 



d q G{z) 



K n 



2=0 

d q logG{z 



dzi 



(16) 
(17) 



z=0 



The counting distribution P(n) is likewise determined by G{z) 

1 d n G{z) 



Pin) 



n\ dz" 



18) 



2 = -l 



Let us now introduce the multidimensional S-variate generating function 



oo oo 



G s ( Zl ,...,z s ) = Yl E •-• E (l + zi) ni ---(l + z s ) ns Ps(ni,...,n s ), (19) 

ni=0ri2=0 ng=0 



from which the S'-variate factorial moments are easily obtained by differentiation: 

/ d Y 1 ( d V s 



H---QS 



n [ ? l] . . . 4 9sl 



dzi 



dz s . 



Gs{z\i ■ ■ ■ , z s ) 



(20) 



z\=— zs=0 



Likewise, S'-variate factorial cumulants are obtained by differentiation of log Gs(zi, . . . , zs)- 

Returning to the function G(z) (§), it is not difficult to see that it can be written 
in terms of the multivariate generating function (|1^) as 

G(z) = G s (z 1 ,...,z s )\ Zi=Z2= ... =zs=z . (21) 

Equation ( |2"T| ) therefore allows to express the factorial moments of n in terms of the 
multivariate factorial moments of {rii, . . . , ns}- 
Application of the Leibnitz rule 



'd_y 

dz 



m = e 



^ ai \a 2 \...a s \ \dz / 



fi(z) 



'd_\ as 
dz 



fs(z) 
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to the function 

f(z) = f 1 (z)...f s (z) 



and using (21) leads immediately to the relation 



The summation is over all sets {a,j} of non-negative integers such that 

s 

E a J = i- 

Formula fl22"|) is a generalization for factorial moments of the usual multinomial theorem. 



F 2 


- 

— r 02 






h 


- -^03 


+ 3(F 1 (2) + F 2 ( 1 2) ) 


_i_ p(2) 
+ -^30 > 


h 


£.(2) 
— -^04 


+ 6F 2 (2) + 4(F 1 (2) 


+ -^31 J + -^40 



Likewise, taking the natural logarithm of both sides of (|2lD , one obtains an identical 
relation as (^) among single-variate and multivariate factorial cumulants. 
As an example, for two groups (S = 2) one finds from (p2|) 

(23) 
(24) 
(25) 

The factorial moments Fq{, Fm, are determined from the counting distribution in a 
single group. The univariate factorial moments F q are obtained from the sum of counts 
in the two groups. The "mixed" factorial moments JPA- (i, j ^ 0) express inter-group 
stochastic dependences. The relations ([2"3"D-([2li)|) are trivially extended to more than 
two groups. 

3 correlations 

3.1 Integral moments and cumulants 

We now apply the general results from the previous section to the case of interest: the 
production of W + W~ in e + e~ annihilation, and their subsequent decay to four jets: 
W + W~ — >• 4g — > hadrons. Here, the number of "sources" S is equal to two. 



Eq. (23) is of particular interest. In a less formal way, it can be written as: 

F 2 = (n{n - 1)) = (ni(m - 1)) + (n 2 {n 2 - 1)) + 2 (nm 2 ) , (26) 

where n\ and n 2 are the number of particles from the decay of W + and W~, respectively. 
This equation could also be derived directly by noting that n = n\ + n 2 and working 
out the expression for (n(n — 1)). To derive relations for S > 2 or between factorial 
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moments of higher order, it is evidently less cumbersome to make use of the generating 
functions and Eq.(^). 

In absence of inter- W correlations of whatever origin, kinematical, dynamical, due 
to experimental selections and cuts, . . . , one has 

( ni n 2 ) = (m) (n 2 ) . (27) 
Exhibiting explicitly the presence of inter-W correlations we write (p6|) as 

F 2 = (n(n - 1)) = (m(m - 1)) + (n 2 (n 2 - 1)) + 2 (m> (n 2 ) (1 + 5/), (28) 

with 

5j = (mn 2 ) I (m) (n 2 ) -1^0, (29) 

a measure of positive or negative inter-iy correlations. Similarly, the factorial cumulant 
K 2 can be written as 

K 2 = K 20 + K 02 + 2 (m) (n 2 > 5 7 . (30) 

This equation expresses the correlation function, integrated over full phase space, of 
the whole system in terms of the integrated correlation functions of each component 
separately, and of the integrated inter-W^ correlation. 



3.1.1 Interlude 

The quantity 5i is related to the variances D w ± of the single- TV 1 * 1 , W ± — ► qq, and 
W + W~ — > 4q multiplicity distributions via the relation 

D ww = ((n - (n})") = ([(m - (m» + (n 2 - (n 2 ))] 2 ) (31) 
= + ^w- + 2 K™2> - 2 (m) (n 2 ) (32) 



giving 

2 (m) (n 2 ) 



n 2 - n 2 - n 2 

8j = ww ^ . (33) 



3.2 Differential distributions 

In section [3.1| , general relations were obtained relating factorial moments and cumu- 
lants of the multiplicity distributions of a W + W~ system to those of its W + and 
W~ components. Since g-th order factorial moments are integrals over phase space 
of g-particle inclusive densities, we now turn to relations among the fully differential 
particle densities and correlation functions of a W + W~ system and those of its com- 
posing parts, considering explicitly possible statistical dependences. We restrict the 
discussion to second-order densities and correlations. 
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For two stochastically independent systems, we derived in an earlier paper || the 
relations 

C 2 WW (1,2) = C 2 W+ (1,2)+C 2 W -(1,2), (34) 
p^ w (l,2) = pf + (1,2) +pT (1,2) + P T(l) P r (2) +p? + (2)p?-(l), (35) 
and further 

p ww (1) = p w +(1) + p w- (1) , (36) 

Here C 2 iVW (l, 2) and (1, 2) are the two-particle correlation functions for W + W~ — > 
4g events and W ± — ► 2g events, respectively; p^ w (l), p^(l) are the corresponding 
single-particle inclusive densities. 

Inspection of ( p0|) suggest to write a general expression for C 2 vVW (l, 2) as 

C 2 ffw (l,2) = C 2 w+ (l,2)+Cr(l,2) + 5,(1,2) {pf + (l) P r(2)+pr + (2)pr(l)}- 

(37) 

The function 5/(1,2) describes correlations among different W*s. 

For 5/(1,2) = (independent W's), (P7|) expresses the additivityf] of the facto- 
rial cumulant correlation functions, a necessary and sufficient condition for stochastic 
independence ||. 

The factors p^ v+ (l)p] v (2), p^ + (2)p^ (1) are introduced for normalization and 
explicitly account for differences in the single-particle densities of the two If' s, as is 
the case when the momentum spectra of the decay products from different Ws are not 
identical i.e. do not fully overlap^. 

With the definitions 

p^ w (l,2) = C 2 ww (l,2)+pf w (l)pf w (2) (38) 
pj(l,2) = C 2 w (l,2)+p7(l)p7(2) (39) 

and using (^TJ) we can write a general form for pY w {l, 2) 

prus) = c 2 w+ (i,2)+cr(i,2)+{pr(i)pr(2)+pr(2)pr(i)}5,(i,2) 
+ {pr + (i)+pr(i)}K + (2)+pr(2)}. («) 

Note that, in general, and in fully differential form, the terms p] v+ (l)p] v (2) and 
p^ + (2) pf (1) are not equal. p^ w (l, 2), however, must be symmetric in its arguments 

3 For independent "sources" , additivity is valid for all orders of the cumulant correlation functions. 

4 In § it was implicitly assumed that pf + (1) pf~ (2) = + (2) ~ (1) = pf(l)pf(2), where 
pY(l) is the single-particle inclusive density of one W and with the further assumption that pY + (1) = 
pY (1) = Pi ? (l). Ignoring possible charge-dependence, these equations are valid if the W + and W~ 
hadronic decay products overlap completely in momentum space. 
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for identical particles. Equation (f4"0|), together with fl39|) takes the form 
pT(l,2) = ^(1,2) + ^"(1,2) 

+ {1+5,(1,2)} { P r(i)pr(2)+p7 + (2)pr(i)} («) 

Defining the experimentally often studied normalized two-particle density 

ww /i n) 

2 [ ' >~ pT(l)pY w (2) , 1 ] 

one finds 

P r(i i 2)+pr(i,2)+{ P r(i)pr(2)+pr + (2)pr(i)} {i+M^g} 
P r (i) P r (2) + P f - (i) P r ( 2 ) + P r (i) P r (2) + P r ( 2 ) P r w 

(43) 



i?^ w (l,2) 



3.3 Correlations in the variable Q 

The previous sections dealt exclusively with fully differential quantities. In practice, 
these are impossible to measure and a projection on a lower-dimensional space is 
needed. We here consider, for illustration, the kinematical variable Q 2 = — (pi — P2) 2 ■, 
the negative square of the difference in 4- momenta of particles 1 and 2. We use the 
notation P w+ (g) P w (Q) for integrals of the type 

/ Jd i p 1 d 3 p 2 p w+ (l)p w -(2)5(Q 2 + (p 1 -p 2 ) 2 ). (44) 

In practical applications, such integrals are calculated using an event and track-mixing 
technique. 

We assume from now on that 5/(1,2) in ( p7|) is a function of Q only: 5i(Q). This 
simplifies the calculations but may not be fully realistic. At least for Bose-Einstein 
studies, it is known that the correlation function of like-sign pairs is not isotropic in 
four-momentum space. To simplify further we also assume that 

pT + ® pT\Q) = pT ® pT{Q) = pT ® p?(Q). 

The expressions in the previous sections take the following form 

P7 W (Q) = p7 + (Q)+p7~(Q)+2{1 + 5 i (Q)} pF®pT(Q), (45) 
C^{Q) = cT{Q) + C^-{Q) + 25 I {Q) {pT ®pTiQ)}- (46) 
Integrating these expressions over all Q, we have 

S ' = (nJ(n w -) J dQS -M{rr®Pr(Q)}- (47) 
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Further 

RWW(0) = P2 + (Q) + P f-(Q) + 2{i + Mg)} {pY + g (g)} 

2 2{pr®^(Q)+Pr®pT(Q)} 
Introducing the "overlap function" 

we can rewrite (^5|)-(|46|) as 

p7 W (Q) = pT{Q)+pT{Q)+2{l + 8 I {Q)}g{Q)pY® P T{Q), (50) 

Cf w (Q) = C¥ + (Q)+C?-(Q)+2Si(Q)9(Q)p7®p7(Q). (51) 

Likewise 

pW w^ _ p7 + (Q) + pf~(Q) + 2 {i + ft(Q)} g(g) g pT(Q) 

l Rf + (Q) , i gQg) , u , , , nU g(Q) 
= 2lT^ + 2i^^ + {1 + W)} iT^g)- (53) 

Here R^(Q) is the normalized two-particle density for a single 

* w » ) = ;4ft» = 1+ ™ (54) 

with K™(Q) the normalized two-particle cumulant density. 

The meaning of the various terms in ([50]) is the following: a term such as (Q) 
"counts" the number of like-sign pairs within a single W + . Integrated over all Q, it 
equals (n w +(n w + — 1)); the term 2 {1 + 5i(Q)} g(Q)Pi ® Pi (Q), "counts" the number 
of pairs with particle i belonging to W + and particle j belonging to W~ . Its 
integral over all Q is equal to 2 (n w +n w -); the term 2g(Q)pf ® pY{Q) = P\ + ® 
Pi (Q) "counts" the number of pairs in uncorrelated Ws. Its integral equals 
2 (^w+) ( n w-)- These relations can serve as a check of any experimental method used 
to calculate the so-called "mixing" terms pf ® pY{Q), pf + ® Pi (Q)- 

To simplify further the expressions, but without essential loss of generality, we shall 
from now on assume that 

pY + (Q) = pY~(Q) = pY(Q). (55) 



We obtain from (53) 



l + #(<?) 1 + #(<?) l + 

i + g{Q) i + g{Q) 
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where K^ W (Q) = i?^ w (<5) — 1, is the normalized cumulant density of the WW system. 

Both K^ W (Q) and K™{Q) are often parameterized with Gaussians. However, from 
Eqs. (|5*B|)-(|57D, it is seen that neither RY W (Q) nor K^ W {Q) will, in general, have the 
same functional Q-dependence as K™(Q) unless K™(Q) ~ 5i(Q) and g(Q) is constant. 

Consider next the following limiting forms of the previous expressions. 

• Fully overlapping and uncorrelated W-decays: <5j(Q) = 0. 

Here, all factors pi<S>pi(Q) are equal so that g(Q) = 1. Equations (|56l-57|) become 

R? W (Q) = l{i + Rj(Q)}, (58) 

K? W (Q) = \k?{Q). (59) 

For this special case, the cumulant correlation function for W + W~ — > 4q events 
is only half that of a single W, corresponding to W + W~ — > 2q events. This result 
was first derived in [§]. Its validity basically rests on the additivity of factorial 
cumulants for sums of independent random variables^. 

• Decreasing overlap, arbitrary Si(Q). Above production threshold, and with 
increasing center-of-mass energy, y/s, the decay products of the two W's will show 
a decreasing overlap in momentum space. The factor p^ + <S> pY {Q) which is a 
measure of this overlap, will tend to zero for any fixed value of Q when s — > oo. 
For g(Q) — > one obtains 

* rw) - a^Sfep™ (61) 

KT"(Q) - K?{Q). (62) 

This result is intuitively clear. The diminishing overlap of the W + and W~ decay 
products decreases the contribution from pairs of particles which originate from 
different W's. For any fixed Q, RY W {Q) increasingly receives contributions from 
pairs belonging to the same W only. 



5 In general, for S fully overlapping identical and independent "sources" it follows from the additivity 
of cumulants that 



_ k^ 



(60) 



(n) q Si- 1 ' 

with kq S ^ the unnormalized integrated factorial cumulant of a system composed of <S* sources, (n) 
its average multiplicity and K^p the integrated factorial cumulant of a single source pj. We note 
that, besides the assumptions mentioned, no further approximations are involved in deriving ( |6C| ) in 
contrast to the derivation in []12| where (pQ) is only obtained under (unnecessary) further conditions. 
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Equations (|56| - |5TD show that, as a consequence of not fully overlapping "sources", 
genuine inter-source correlations, represented by Si(Q), are "diluted" in the actual 
measurement of R™ W (Q) or (Q) since g(Q)/(l + g(Q)) may become quite small 
in practical situations. 

This has further implications. In experimental conditions where one is attempting 
to avoid the effect of inter-source correlations (such as in VT-mass measurements), 
methods should be devised to render the overlap function as small as possible e.g. by 
suitably chosen cuts. 

On the opposite, when the main interest is focussed on establishing inter-source 
correlations, one could try to devise methods that maximize g(Q)/(l + g(Q)) and thus 
increase the experimental sensitivity to Si(Q). 

With respect to searches for inter-jy correlations at LEP, it remains to be inves- 
tigated how much the presently adopted cuts and algorithms used to select with high 
efficiency WW events, affect the shape of the overlap fuction and contribute to weaken 
(possible) genuine WW correlation effects in the analyzed data. The same holds for 
studies of possible color-reconnection effects. 

With the formalism described above, it is now straightforward to obtain expressions 
for a variety of observables used in experimental studies of inter- W correlations. Some 
of these are reviewed in the following. Others are easily constructed using the results 
of this section. 



4 Examples 

4.1 The L3 test-statistics D and D' 

The L3 collaboration, in a search for inter-11 7 (BE) correlations, discussed in a recent 
paper [[13] the distributions ("test-statistics") 



Ap(Q) = pf w (Q)-2pJ(Q)-2 P r x (Q), (63) 
2pY(Q) + 2 P r*(Q) 



D(Q) = H2 , (64) 



D\Q) = D(Q) ; (65) 

^MC, no BEC(V) 

where P2 11X (Q) is identical to pf + ® pf {Q)- The ratio D'(Q) is obtained by divid- 
ing D(Q) by the same derived from a Monte Carlo calculation without Bose-Einstein 



correlations. As very similar quantity was studied by ALEPH [14J . 
These distributions can be rewritten as 

Ap(Q) = 2g(Q)5 I (Q)pY®pY(Q), (66) 
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D[Q) = 1+MQ) wwhry (67) 



Observe that (S7j) and ( 68f) are formally, but not necessarily experimentally, identical^. 

These expressions show that none of the studied observables isolates completely the 
genuine inter- W correlation function 5i(Q). It can be measured most directly via the 
unnormalized function Ap(Q), as seen from (|66"D, once the mixing term ® (Q) 
is determined. This equation therefore (strongly) suggests to study the ratio 

A " (Q) -*,(«). (69) 



2 P r x (Q) 

It can be determined directly from data only 



4.2 The fraction of pairs from different W's 

This distribution was used by L3 |L5| and DELPHI [|l^]. The fraction of like sign pairs 
from different W's is, in our notation, defined as 



2{1 + 6 i ( Q)}pY + ®pY~(Q) 

pI w (Q) 

or 



j?(n\ ^l"- ^ u iW)s Pi w Pi / 7n x 
*[Q) = > l 7U J 



F(n) = g(Q) {1 + 5j(Q)} 
W) Rf(Q)+g(Q){i + Si(Q)Y 1 } 



The equation can be solved for g(Q) with the result 



^'^ i-lw^w (72) 



Using -F(Q) presented in fj"5| , we have plotted and g(Q)/(l + g(Q)) in Fig. 1. 

Since 5/(<5) is unknown, we assume, for illustration, the "Goldhaber"-form Si(Q) = 
Aexp(-rfQ 2 ) and Rf = 1 + Aexp(-r 2 Q 2 ), with r = 0.67 fm and A = 0.7 taken 
from [y|. We set rj = r. 

The function g(Q) /(l + g(Q)) has a value between 0.14 — 0.24 at Q = 0, depending 
on A, implying (c/r. Eq. (|57D ) that the effectively measured strength of possible inter- 
im correlation effects is smaller by the same factor. 

With g(Q) derived from the L3 results, we calculated the quantity D(Q) as a 
function of A. The result is shown in Fig. 2. As is clear from (|67D , D(Q) is not of the 
form ~ (1 + Aexp (-k 2 Q 2 )), used in jEJ. 

6 Eq.(|6^) remains unchanged when it is calculated from a Monte Carlo with only intra- W correla- 
tions. 
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To test for inter- W correlations, a density function, R^ w (Q) B es: is defined in |T5| 
which is supposed to coincide with the data distribution R^^(Q) (defined below) 
if inter- W correlations are absent. This function is stated to describe the data on 
RYl3(Q), thus adding extra support to the claimed absence of inter- W BEC. 

The quantity RYl3(Q) is defined as the ratio of the like-sign two-particle density 
in the data to the same obtained from a Monte Carlo model without Bose-Einstein 
correlations. Assuming that the model reproduces perfectly RY(Q), this quantity is 
given by 

R ww fn , RJ + g(Q) {1 + 5j(Q)} K^ + g(Q)6 I (Q) 

The "test" function R^ W (Q) bes is taken to be 

< W (Q) BES = 1 + {1 - F(Q)} K?(Q), (74) 
with F(Q) given by (|7ll) and setting Si(Q) = 0. Equation ([74]) can then be written as 

^-''^ i^Rqi ^ (75) 



With 5j(Q) = in fl73|) for no inter- correlations, it is seen that this distribution 
and that of Eq. ( |75| ) are different, unless K™(Q) = or g(Q) = 0. We conclude that 
(|T1D is not suitable to test for the "null hypothesis" Sj(Q) = 0. An analogous test- 
distribution, used by the DELPHI Collaboration for similar purposes (see below), has 
the same defect. 



4.3 The DELPHI observables 



In the DELPHI analysis of Bose-Einstein correlations in WW events |16|], the following 
ratio is studied ww 

D data _ P2 

^2,MC, no EEV*) 

where the denominator is calculated from a Monte Carlo without Bose-Einstein effects. 
This can be rewritten using ( (50|) as 



which coincides with the general expression fl56|). We have assumed that, in the Monte 
Carlo, identical particle pairs are uncorrelated, i.e. pY {Q) — pT ® pY(Q)- 

The function ( |77| ) is plotted in Fig. together with R^ (denoted by i?2q in [0) 
measured in the DELPHI analysis and used here as input. From the above expression, 
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and from the figure, it is clear that this distribution remains sensitive to inter-H^ 
correlations. However, it cannot, in general, be described by the same parameterization 
as used for R™(Q), even for 5j(Q) = 0. Fig.^ also shows, for illustration, the same 
function using a Gaussian parameterization of $i(Q) with a "radius-parameter" of 
3 fm. The latter might be typical for the range of second-order interference (HBT) 
correlations associated with incoherent strings. In this case, the Q-region affected by 
BEC is restricted to values below ~ 100 MeV and will decrease even further for still 
larger "radii". 

In Fig. | we show the intercept, R^ ta (Q) at Q = 0, and that of R^(Q), as a 
function of A. It is seen that the intercept for the WW — > 4q channel remains below 
that of RY(Q) for A < A. It is, of course, independent of the values of the parameters 
r and 77. 



4.3.1 i?4q(Q) (mixing) 

DELPHI also studied a "test" distribution, R^qiQ) (mixing), constructed from mixed 
independent (WW — ► 2q) events and defined as 



R 4q (Q) (mixing) = r — l - Jdata (78 







(Q) 


data 


pf(Q) + P T*{Q) 


MoCa, No Be 



with p™ lx (Q) = pf + ® pf ■ The above expression is equal to 

i? 4q (Q)(mixing) = l + T ^i^. (79) 

This expression indeed coincides with ( fTTD for Si(Q) = and can thus be used to test 
for inter- W correlations. It is plotted in Fig. |3] and depends here (weakly) on Sj(Q) 
which enters the definition of g(Q). We remark again that R 4( ^{Q) (mixing) cannot, in 
general, be described by the same parameterization as used for R^(Q). 

Defining AA(mixing) as the difference in intercepts at Q = of R^q(Q) (mixing) and 
^data we find 

AA(mixing) = A 9 ^°] , (80) 
V S; 1 + 0(0)' 



with A = 5/(0). 



4.3.2 i*4 q (Q) (linear) 



A further quantity studied by DELPHI is R^qiQ) (linear), defined in Eq. (25) of ref. |L6|. 
It is constructed to supposedly coincide with R^ W (Q) if inter- W correlations are ab- 
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sent. With our notation it is given by 



i? 4q (Q) (linear) = itf (Q) + {l - Rf(Q)} ^^J— . (81) 

This function is plotted in Fig. |5|. The distribution differs from R^" 3 " even for 5i(Q) = 
(or A = 0) and, therefore, has not the intended properties. Nevertheless, from Fig. 3 
and Fig. 5 it can be seen that R^iQ) (linear) and i?4q(Q) (mixing) are numerically quite 
close. The weak dependence on A is again due to the apriori unknown dependence of 
g(Q) on 5j(Q). 



5 Summary 

In this paper we described a mathematical formalism which should allow a systematic 
study and improved understanding of particle correlations in a physical system which 
is composed of S possibly stochastically correlated parts. Most emphasis was put on 
two-particle inclusive densities and correlation functions for S = 2, but the formalism 
can be extended to arbitrary orders and arbitrary S. 

Taking the reaction e + e~ — > W + W~ — > qi^QzQi ~^ hadrons as a prime example, 
it was shown how to relate the measurement of correlations within one hadronically 
decaying W to those measured in the full WW — > 4q system. 

For the observables used in present studies, correlations in fully hadronic WW 
decay's can be expressed in terms of the correlation function of a single W, an inter-W 
correlation 5i and a function, g, which quantifies the degree of overlap in momentum 
space of the decay products from different Ws. 

On several examples, we have shown that some of the presently used techniques to 
search for inter-H^ correlations are not optimal. In some cases, experimental quanti- 
ties have been used which lack the correct mathematical properties and are based on 
intuition rather than on sound mathematics. 

The determination of the overlap function allows to assess quantitatively the sen- 
sitivity of a particular distribution to inter-1^ correlations. Moreover, the influence of 
experimental and methodological cuts can be systematically investigated . 

An estimate of g{Q), presented here, indicates that the fully hadronic, four-jet decay 
channel of the reaction e + e~ — > W + W~ may not be optimal to search for inter-VT Bose- 
Einstein effects. The usual event selections, requiring four well-isolated jets in the final 
state, weaken the sensitivity to possible inter-W^ correlations. Better results may well 
be obtained in analyses using less stringent cuts, or which include e.g. three-jet events 
with suitably chosen topology so as to keep the QCD background at a tolerable level. 

Our survey of some of the test-distributions used so far indicates that the observable 
Ap(Q), defined in Eq. (|63"D, proposed in and measured by L3 fT3]| , provides a direct 
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measure of possible inter-VF correlations. It suggests itself as suitable observable well- 
adapted for a global LEP-wide combination of WW data. 
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Figure 2: The test-statistic D{Q) used by L3, Eq. fl67|) . Also shown is the L3 parame- 
terization of R2q(Q) = RY(Q)- 
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Figure 3: DELPHI: The normalized densities R 4q {Q) = R^ W (Q), Eq.(^) in the 
text, using as input the DELPHI parameterization of i?2 3 (Q) = RY{Q) an d assuming 
$i(Q) — Aexp (— r 2 Q 2 ); A varies in the range 0.0 — 1.0, in steps of 0.2 The dashed 
curves show the same with r = 3 fm. Also shown (dot-dashed) is the distribution, 
Eq.(|79|)) corresponding to mixed independent WW — > 2q events. It depends implicitly 
on Si(Q) through the function g(Q). 
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Figure 4: DELPHI: The value of the normalized density R 4q (Q) = Rj w (Q), Eq.([F7|) 
in the text, and of R2 q (Q) = R^iQ), at Q = 0, as a function of the parameter A in 
6j(Q) = Aexp(-r 2 Q 2 ). 
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Figure 5: DELPHI: The normalized density R Aq (Q) = R^ W (Q), Eq.(^) in the text, 
using as input the DELPHI parameterization of -R 2 q(<3) = R™(Q) and with 5j(Q) = 
A exp (— r 2 Q 2 ); A varies in the range 0.0 — 1.0 in steps of 0.2. Also shown is the 
distribution corresponding to the so-called Linear Scenario in ref. [16 
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